Soliton appearing in boson-fermion mixture at the third order of the interaction radius 
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In this paper we consider an ultra-cold mixture of boson and fermion atoms on the basis of quantum 
hydrodynamics. Small perturbations in such systems are being analyzed. A possibility is shown for soliton 
solutions of a new type to appear if the third order of the interaction radius is taken into account in the 
analysis of interactions. A fermion-fermion interaction occurs in explicit form if this approximation is 
accepted. The conditions that lead to occurrence of this type of soliton in a mixture of boson and fermion 
atoms were investigated. Restrictions on the fermion-fermion interaction were found that are necessary 
for this kind of perturbations to appear in the system. Conditions determining whether perturbances would 
be a condensed soliton or a rarefied soliton are shown. Requirements to the experimental detection of a 
new soliton type in boson-fermion mixture are considered. 
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I. INTRODUCTION 

Experimental and theoretical approaches to the physics 
of quantum gases have been continued to grow in recent 
times, and various nonlinear structures has been studied 
there |l]]- 01 ■ Methods for producing Bose-Einstein con- 
densate and investigating its properties have been devel- 
oped (SO. This showed stable conditions to exist both for 
repulsion and attraction between the atoms, for example, in 
a Bose condensate of 7 Li atoms J^]. A mixture of Bose 
and Fermi degenerate gases, where interaction between 
atoms may be tuned from repulsion to attraction, is consid- 
ered in nXOn - Papers IIIOL fllll have been devoted to investi- 
gations of ultra-cold boson-fermion mixtures. A branch of 
physics that works on producing of a stable boson-fermion 
mixture is being developed as well lfl2ll . 

A lot of interest has been attracted in recent years to 
the investigation of soliton solutions. Papers devoted to 
properties of different types of solitons are constantly be- 
ing emerged. As to the present time a bright soliton has 
been produced in a Bose-Einstein condensate with a re- 
pulsive interaction between atoms. A soliton of this type, 
has been produced, for example, in the condensate of &7 Rb 
atoms in a weak periodic potential lfl3ll . The formation of 
a bright soliton in the condensate of 7 Li atoms is consid- 
ered in J^HH)] for a quasi-lD case. Properties of a dark 
soliton have been studied in 111 511 . They show experimental 
generation of a dark soliton in a cigar-shaped condensate 
of 87 Rb atoms with the use of phase imprinting technique. 
The study of a dark soliton in Bose-Einstein condensate 
inside a vortex ring is described in iflrjll . In addition to 
one-dimension solitonic perturbations in a Bose-Einstein 
condensate, three-dimension perturbations are also attract- 
ing attention lfl7ll . A possibility was shown for solitonic 
perturbations to exist in a Bose-Einstein condensate with 
dipole-dipole interaction lfl8ll . In recent years predictions 
have been made that solitons of bright-bright, dark-bright 
and dark-dark types may exist in a two-component BEC 
ifl^l . Solitons in boson-fermion mixture is special and in- 
teresting item. Perturbations in boson-fermion mixtures 
have also been studied and the existence of a bright soli- 
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ton in boson-fermion mixture has been proved in [10, 20]. 

In the present work we study ultra-cold boson-fermion 
Examples of such mixtures are 7 Li — 6 Li 
6 Li JH and 87 Rb - 40 K [24]. The pos- 
sibility of new solitonic perturbations in a mixture of bo- 
son and fermion atoms with short-range interaction po- 
tentials is analyzed in particular. Theoretical analysis at 
higher order of interaction account accuracy allows predic- 
tion of new types of soliton solutions 112511 - This approxima- 
tion can be derived with the aid of quantum hydrodynam- 
ics approach, which has been developed in recent decade 
1 25142711 . The set of quantum hydrodynamic equations de- 
scribing an ultra-cold mixture of boson and fermion atoms 
may be derived from a many-particle Schrodinger equa- 
tion ll25ll . In |28ll . for example, they used this approach to 
study how the shape of a bright soliton changes as inter- 
actions in a Bose-Einstein condensate are considered more 
precisely. The method of quantum hydrodynamics is very 
useful in different areas of physics. Namely, it was quan- 
tum plasma ll2r3l . J2^], plas ma of particles with the own 
mag netic moment IDOfl - 13411 . relativistic quantum plasma 
113511 . ultracold Bose and Fermi gases with nonlocal inter- 
action J25I], quantum particles with electrical ll36ll . ll27ll 
and magnetical 13711 p olarization, particularly being in the 
BEC state j27il, ! 137h, the graphene electrons S and 
BEC of graphene excitons II39II . 

Most of theoretical papers devoted to the ultracold 
fermions contain nonlinear Schrodinger equation (an ana- 
log of Gross-Pitaevskii equation using for bosons) which 
does not include interaction ll40ll - ll49ll . It is connected with 
the fact that the first Born approximation gives no contribu- 
tion in interaction of ultracold fermions. This equation un- 
doubtedly accounts the Fermi pressure contribution caused 
by Pauli principle. It is very interesting to understand a 
role of interaction in fermion systems. For example, vor- 
tices in a strongly interacting gas of fermionic atoms on the 
BEC- and the BCS-side of the Feshbach resonance was ex- 
perimentally studied in Ref. [50]. Presented in our paper 
model based on the quantum hydrodynamics method con- 
tain fermion-fermion interaction whose form was derived 
directly from the Schrodinger equation QZSJ. 
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In this paper we discuss the existence of new types of 
solitons in a mixture of boson and fermion atoms and deter- 
mine the conditions of their experimental detection. A sub- 
stantial role in our analysis is given to the fermion-fermion 
interaction coefficient, which appears in the analysis of the 
third order of the interaction range. A perturbation method 
introduced by Washimi 115 ill . 115 211 is used here to find soli- 
tonic perturbations. An analogous method was used for 
a spinor-1 BEC studying 115 311 . Other methods for weak- 
nonlinear analysis of BEC have been considered in litera- 
ture. For example, in Ref. [54] the Krylow-Bogoliubov- 
Mitropolskii method was used for nonlinear frequency shift 
calculation. It was shown that account of interaction up 
to TOIR approximation leads to the new solitons in BEC 

Our paper is organized as follows. In Sect. 2 we present 
basic equation and describe using model. In Sect. 3 we 
study solitons in boson-fermion mixture which appearing 
due to interaction account up to TOIR approximation. In 
Sect. 3 detailed analysis of conditions of a "fermion" soli- 
ton existence is presented. In Sect. 4 condition of a "bo- 
son" soliton existence is described. In Sect. 5 brief sum- 
mary of obtained results is presented. 



II. THE MODEL 

Let's consider a mixture of ultra-cold boson and fermion 
atoms having our focus on one-dimension perturbations. 
The set of quantum hydrodynamic equations consists of the 
following equations: the continuity equation for bosons 

dn b d(n b v b ) 
dt ox 
the momentum balance equation for bosons 

dv b 1 dvi 2 h 2 d ( 1 d 2 Jn~ b 

m b n b — + -m b n b — — n b — —— v 2 

at 2 dx 2m b ox \ ^JrT b dx 



dn h 



d 3 



i>7. 



ry. 

- J- bbn b ~^ — 1 2bb - 3 

Ox 16 dx 
„ dn f 1 d 3 n f 



the continuity equation for fermions 

dn f d(n f v f ) _ 



+ 



0, 



dt dx 
and the momentum balance equation for fermions 



(3) 



dv f 1 dvl 
m f n f —r L -\ — mtnr— 1 - 
1 1 dt 2 1 1 dx 2m f 1 dx \ ^/nj dx 



h 2 n d_( 1 d 2 ^nj 



3 dn f dn 2 3 



8 2// dx dx 8 



2//' 



d 3 n f 
'dx 3 ' 



1 dn* 8 / 3 
+ 2 (3?r) T2ff ^x~ + 



5m 



■(3tt : 



2x2/3 dn f 5 / 3 



dx 



~ dn b 1 d 3 n b 
Tbfnf ^x- + 2 T ^ n ^' 



(4) 



The equations given above utilize following designations: 
m b , nif- masses of boson and fermion atoms, respectively; 
n b , rif - concentrations of bosons and fermions; v b , Vf - 
respective velocity fields. T bb , T b f - coefficients of boson- 
boson and boson-fermion interactions at the first order of 
the interaction range. T 2bb , Y 2b f, T 2 // - coefficients of 
boson-boson, boson-fermion and fermion-fermion interac- 
tions at the third order of the interaction range. Interaction 
coefficients can be defined with following equations: 



fay — — 



4tt 

T 

4tt 
15 



drr 3 



dU l3 (r) 
dr 



dr 



(5) 



(6) 



Presented here hydrodynamics equations correspond 
to the system of two nonlinear Schrodinger equations 
one for bosons another for fermions [|25J. These non- 
linear Schrodinger equations are nonlocal and integro- 
differential. The nonlinear Schrodinger equation for Bose 
particles is a generalization of Gross-Pitaevskii equation. 
This generalization appear due to more detailed account 
of short-range interaction up to third order on interaction 
radius. For the first time a nonlocal Gross-Pitaevskii equa- 
tion were derived in Ref.s. [|56ll . n57n . This equation does 
not contain integral terms. Brief comparison for different 
nonlocal generalization of the Gross-Pitaevskii equation is 
discussed in ref. Il28ll . 



Coefficients T bb and T b f are related to interactions co- 



efficients 
follows: 



and g b f of the Gross-Pitaevskii equation as 



T 



bb 



T 



bf 



-9b j ■ 



(7) 



It should be noted that the following relationship exists 
between interaction coefficients and scattering amplitude 

4irh 2 an 

Ty = (8) 

m 

All equations above are approximated up to the third 
order of the interaction range. If terms containing T 2 b&, 
^ ' ~^2&/> ^2ff > are neglected then the resulting equation set 
corresponds to the first-order approximation (the Gross- 
Pitaevskii approximation) which has been used, for exam- 
ple, in 111 ill . It was shown in Il25ll that the approximation 
is possible that derives coefficients T 2 bb, ^2bf from T bb , 
T b f. The explicit form of this approximation is: 



T ohh ~ r- 2 T 



and 



2y 



6/: 



where r - is a constant of the same order of magnitude as 
the atomic radius. 
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HI. THE SOLITON SOLUTION IN A BOSON-FERMION 
MIXTURE. 

The equation set in question should be expected to have 
soliton solutions of a new type due to more accurate recog- 
nition of atomic interactions. The method of perturbations 
may be applied to find this soliton I15lll58ll . According to 
this method all hydrodynamic values may be represented 
as: 



n b = n ob + enib + e n 2b + ••• , 



n } = n of + en lf + e n 2f + 



v b = ev lb + e v 2b + 



(9) 
(10) 



(11) 



Vf = ev lf +s 2 v 2f + .... (12) 
We also performed the following "scaling" of variables: 

£ = e 1/2 (x-Ut) (13) 

and 

T = e 3/2 Ut. (14) 

The latter expression introduces so-called "slow" time. 

The following expression of the phase velocity U can be 
derived from equations (fl]|4| in the first order of the small 
parameter e: 

— 777 jTlpkT bb ) 



ul 



2m b mf 



1 



± 2m m ^ ( 0mfcn o/ ~ rn f n 0b T bb j 2 



+ 4m b m f n 0b n f(®T bb + T 



bf. 



1/2 



(15) 



where 



e 



(3vr 



2\2/3 



K 2 



n ^ \3rrif 

Quantity describes the contribution of fermion dynam- 
ics in boson-fermion mixture. O consists of two parts. The 
first one appears from Fermi pressure and the second term 
presents the fermion-fermion short-range interaction which 
arises at interaction account up to the third order of the in- 
teraction radius. In formula (Q3) the quantity under radical 
is positive. 

Since our mixture is comprised by two interacting sub- 
systems, boson atoms and fermion atoms, two different val- 
ues of phase velocity are obtained here. We call phase 
velocity U + which has "+" sign before the radical, the 
"fermion branch" of the solution. Phase velocity U- is 
called the "boson branch" of the solution by analogy. We 
introduced these terms to emphasize the fact that in the ab- 
sence of boson-fermion interactions T b f = the system 



can be divided in two subsystems of bosons and fermions, 
respectively, which do not interact. Phase velocity U + cor- 
responds in this case to fermions 

Ul = 

m f 

and U- corresponds to bosons 

U 2 -- 



n ob 

1 66- 



m b 

We can see that U 2 > for particles with repulsive inter- 
action T bb < 0. 

In the second order of the small parameter e we derive 
the Korteweg - de Vries equation for small perturbations of 
boson concentration in a boson-fermion mixture: 
3r7 lb d 3 n lb dn lb 

P ^r + q ~de + snib ^r = °' (17) 

where coefficient p of the term containing slow time 

p = 2U 2 (<dm b n 0f - T bb n ob m f - 2U 2 m b m f ), (18) 
coefficient q of the term that corresponds to dispersion 

q = — (@n f — mfU 2 )x 

h 2 1„ 1„ T bb 1„ U 2 m 



~. 1" -T 2 hb77ofc — -'^2bfn'Q b — — 

4m f 8 2 T bf 



^2bf~^ 

1 L bf 



+ (im7 + I T2//ri0/ ) ( n ° fcTbfa + u2m b) 

and coefficient s of the non-linear term 

U 2 m b 



(19) 



n ob 



-(Qnof-U m f ) + T bf (T bb n ob + U m b ) 



U 2 m b 



+ r bb + ^^ ^ 



n ob 



n ob 



T 



bb 



20 



2\2/3 



73 
0/ 



(20) 

The soliton solution of the Korteweg - de Vries equation 
(fTTl) is well known in the form of 

3pV 1 



77 



1ft 



coshZ ( Jiiv 



(21) 



where 7/ = £ — Vr. 

Below we discuss perturbations in bosons implying that 
perturbations in fermions are similar to them due to the 
following linear relationship between n lb and n\f. 



T hf n 



6/™!/ 



T b6 H | n lb . 

n ob 



(22) 



Soliton pairs (in bosons and fermions, respectively) of 
two types may evolve: dark-bright and dark-dark. 
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FIG. 1. (Color online) The relationship between dimensionless 
phase velocity Wf and coefficients of fermion-fermion interac- 
tion y b f and boson-fermion interaction -y b f for fermion branch. 
Figure is made for the following parameters of the system: 7^ = 



-1CT 



1, N = 1, B, 



bf 



io~ 



and B, 



66 



-1CT 



This 



plot and plots 2-4 made for the repulsion between Bose particles 
~f bb < 0. Dimensionless phase velocity Wf has analogous form 
in the case of attraction between Bose particles. 



IV. THE FERMION BRANCH OF THE SOLUTION. 

Let's investigate "fermion branch" of the solution, where 
phase velocity U from the equation (TT3T > is taken with " + " 
sign. 

If boson-fermion interactions are "strong" enough com- 
pared to boson-boson interactions, then U 2 > 0. We use 
apply assumption to the fermion branch below. 

For numerical analysis of our results we introduce fol- 
lowing dimensionless variables: dimensionless phase ve- 

1 /3 

locity W+ = nifU + /(hn Qf ), mass rate fi = mf/m b , 
concentrations rate N = y/n ob /n f, and dimensionless 
interaction constants 7 fcb = mfTi 0b Y bb /(nQf 3 ti) 2 , 7b/ = 
mf^n^jThf/inlfh) 2 , j ff = m / n 0/ T 2// /n 2 , 



B b b = mfn 0b T 2b b/h 2 , and B, 



bf 



n 



2/3 

of 



bf- 



Phase velocity for fermion branch is shown at Fig. 1 . 
It can be shown that 



p < 



(23) 



at any allowable values of physical parameters. It follows 
from numerical analysis of p. 

As a starting point we considered the first order of the 
interaction range, neglecting the contribution of terms that 
are proportional to coefficients T 2bb , ^2tf, Y 2// in equa- 
tions © and ©. 

In this case the coefficient q takes the following value: 



If 



a- 



4mi 



h 2 



Am f 



>0, 



(24) 



|y b f| 




0.00 

-0.01 [1] 
-0.02 * 



FIG. 2. (Color online) The relationship between soliton ampli- 
tude H b of Bose subsystem and the coefficient of boson-fermion 
interaction 7^/ and fermion-fermion interaction jtt for fermion 
branch. Figure is made for the following parameters of the sys- 



tem: 

B bb 



766 = 
= -HT 4 



-1CT 



H = 1, JV = 1, B bf = 10 



-3 



and 



where a > and /3 > 0. Evident form for a and /3 to be 
1 



a 



2m b 



(Qm b n of + T bb m f n ob ) 



+ \l (Qm b n of + T bb m f n 0b ) 2 + 4m b m f n 0b n 0f (T bf ) 2 

(25) 



and 



1 



2m i 



(Qm b n of + T bb m f n 0b ) 



(@m b n of + T bb m f n 0b ) 2 + 4m b m f n 0b n 0f (T bf ) 

(26) 

Expressions (l2"T"ll . (1231 and (l24l make it clear that no 
soliton exists in the boson-fermion mixture in the Gross- 
Pitaevskii approximation, as the radicand in d2J~b is nega- 
tive. 

So, we took into account the third order by the interac- 
tion radius. 

Now the coefficient q has the following form: 



q 



P 4m t 



+ 



h2 1 v 
h -n 0f Y 2f f 

4m f 4 1 11 



+ 



(27) 



It means that if the repulsion of fermions is strong enough 
(i.e. the coefficient of fermion-fermion interaction Y 2 // < 
0) the q value would be negative. Strong fermion-fermion 
interactions in the boson-fermion mixture make possible 
perturbations of a new type, which do not occur in the first 
order of the interaction range. 

In the case the considered perturbations of concentra- 
tion are small, solitons of two kinds may arise: solitons 
of fermion rarefication and of fermion compression. The 
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FIG. 3. (Color online) The dependence curve of the soliton am- 
plitude Ef of Fermi subsystem vs the coefficient of fermion- 
fermion interaction jff and boson-fermion interaction y b f for 
the fermion branch. The plot is built with following parameters: 
766 = -1CT 2 , fi = l,N = l,B b f = 1CT 3 , and B bb = -1(T 4 . 




FIG. 4. (Color online) The dependence curve of the soliton am- 
plitude Ef of Fermi subsystem vs the coefficient of fermion- 
fermion interaction 7/ f and boson-fermion interaction y b f for 
the fermion branch in the case negative y b t that corresponds to 
repulsion between Bose and Fermi particle. The plot is built 
with following parameters: ^ bb = — 10 , fj, = 1, N = 1, 
B bf = 10 -3 , and B bb = -10 -4 . 



soliton type is determined by the sign of its dimensionless 
amplitude of bosons 

3pV 



■=■6 



sn ob 



(28) 



and fermions S * which connected with the H b by formula 



The dependence between amplitude of soliton in the bo- 
son subsystem and the coefficient of boson-fermion inter- 
action T b -f is shown at Fig.s 2-4. If the coefficient value is 
high positive (it corresponds to a strong attraction between 
boson and fermion atoms) then boson perturbations would 
be solitons of depression (dark soliton, it is shown at Fig. 2) 
and, as it follows from (l22l . fermion perturbations would 
be solitons of compression (bright soliton) for j b f > or 
soliton of depression (dark soliton) for 7 b j < (it is shown 
at Fig.s 3 and 4). These are dark-bright and dark-dark soli- 
tons in a two-component system of bosons and fermions. 
So, a strong repulsion between bosons and fermions (i.e. 
negative values of the T b f coefficient) would lead to soli- 
tons of compression of bosons and fermions, which corre- 
spond to a dark-dark soliton. 

The dimensionless soliton width D can be defined by the 
following expression: 



D 



Vp 



(29) 



It follows from (fl"8l and (l27l l that soliton width depends 
on the fermion-fermion interaction. This dependence is 
plotted at Fig. 5. Note that the surface's shape does not 
change along with the sign of j b j as the expression for the 
soliton width D contains the coefficient of boson-fermion 
interaction in a squared form. 

There is interesting behavior of this solution for attrac- 
tive boson-boson interaction. In this case soliton exist 
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-40 _ 30 
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FIG. 5. (Color online) The dependence curve of the soliton width 
D vs the coefficient of fermion-fermion interactions 7// and 
the coefficient of boson-fermion interaction ^ b * for the fermion 
branch expressed in reduced values. The plot is built with follow- 
ing parameters: j bb = — 10~ 2 , fi = 1, N = 1, B b j — 10~ 3 , and 
B bb — — 10~ 4 . Analogous dependence take place at ~f bb > for 
Ibb ~ 10~ 2 -r 10~ 3 in the area of large boson fermion interaction 
| 7 fc/ 1 > 10. 



in two areas of system parameters. For the large \j b f\ 
(|76/l > 10)' function -D(7//,76/) is analogous to the 
same presented at Fig. 5. But for 7^ > there is an area 
at the small {^fl where exist soliton solution. In this area 
soliton width is large, much more than 1. It is presented at 
Fig.s 6-8. 

In the case j bb > and \ j b f\ < 1 amplitude of bosons is 
positive E b > 0. So, we have bright soliton in subsystem 
of bosons Fig. 6. From Fig. 7 we can see that if 75/ < 
amplitude of soliton in Fermi subsystem is positive. Thus, 
we have found bright soliton in Fermi subsystem. If 7^ > 
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FIG. 6. (Color online) The dependence curve of the soliton width 
D the coefficient of boson-fermion jb t and fermion-fermion jff 
interaction for the fermion branch. This plot made for the attrac- 
tion between Bose particles 7 > 0. The plot is built for two values 
of 7f,b (for blue-narrow surface we took j b b — 10~ 3 , and for red- 
wide surface we used 7^^ = 10 -2 ) with following parameters: 
p = l,N = l,B b f = 1(T 4 , and B bb = -1CT 4 . 

we have dark soliton in Fermi subsystem. In the result 
in boson-fermion mixture we have bright-dark soliton for 
7b/ < 0, and bright-bright soliton for j b f > 0. 

Considered in the paper values of interaction parameters 
may be obtained using Feshbach resonance ll59ll60ll . It fol- 
lows from (|5]l and ©, that changing potential of atomic 
interaction changes both Ty and T 2 ij coefficients. Note 
that soliton of this type occurs due to boson-fermion inter- 
actions. 



V. THE BOSON BRANCH OF THE SOLUTION. 

To obtain the boson branch phase velocity U from the 
equation ( fT3T > should be taken with "— " sign. 

Using the first order of the interaction radius does not 
allow to predict solitonic perturbations in boson-fermion 
mixtures. It can be done only when the third order is taken 
into account. If the mixture of bosons and fermions di- 
vides into two subsystems that do not interact, bosons "do 
not feel" the presence of fermions, T h f = 0, and the coef- 
ficient q ~ — fi 2 /4mf, — n 0b T 2b b/8 115511 . This means that 
repulsion of bosons (T 2bb < 0) must be strong enough to 
cause soliton formation in a subsystem of boson atoms. 

If the interaction of bosons and fermions occurs, then 
the acceptable range of physical values, where the soli- 
ton formation is possible, changes. In boson-fermion mix- 
ture a soliton corresponding to boson branch exists at more 
large strength of boson-boson interaction, in comparison 
with the case when T b f = 0. We have bright soliton 
solution (H;, > 0) in Bose subsystem. In Fermi subsys- 
tem we find dark H/ < (bright 3j > 0) soliton for 
repulsion (attraction) between Bose and Fermi particles 




FIG. 7. (Color online) The dependence curve of the soliton width 
D (blue-top surface) and amplitude for bosons E b (red-middle 
surface) on the coefficient of boson-fermion jbf an d fermion- 
fermion jte interaction for the fermion branch. Soliton exist if 
its width D is positive. We present soliton width D and amplitude 
for bosons on the same plot to show the value of amplitude in 
the area where soliton exist. Green plane-lower surface present 
zero level. It show us that amplitude is positive. The plot is 
built with following parameters: 7 fcb = 10~ 2 , y, = 1, N = 1, 
B bf = 1(T 4 , and B bb = -1(T 4 . 




FIG. 8. (Color online) The dependence curve of the soliton width 
D the coefficient of boson-fermion j b f an d fermion-fermion 7// 
interaction for the fermion branch. The plot is built with follow- 
ing parameters: jbb = 10~ 2 , [i = 1, iV = 1, Bbf = 10~ 4 , and 
B bb = -10~ 4 . 



T b f < (Tbf > 0). We have got it because quantity 
— (T fcb + m b U^_/n b) presented by formula (l22l is posi- 
tive. Consequently, Sf has the same sign as T b f. So, in a 
mixture of bosons and fermions this type of soliton solution 
may exist. 

Finally, we can note that two soliton solutions exists in 
boson-fermion mixture due to account of short range inter- 
action up to the third order of the interaction radius. 
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VI. CONCLUSION 

The quantum hydrodynamics approach is well suitable 
to derive equations that describe multiparticle quantum 
systems. In our case we applied the QHD approach to the 
ultra-cold mixture of boson and fermion atoms with short- 
range interaction potential in order to analyze the possi- 
bility of small-scale solitonic perturbations. We show that 
solitons of two novel types may occur, which are related to 
fermion and boson branches of the solution, respectively. 

In the fermion branch the soliton of a new type occurs 
due to boson-fermion interactions at the first order of the 
interaction range. Fermion-fermion interactions also play 
substantial role in formation of these perturbations. The 
type of soliton formed in the system (compressed or rar- 
efied) is determined by the sign of boson-fermion interac- 
tion, i.e. repulsion or attraction, respectively. The depen- 
dence of fermion concentration on system characteristics 
was successfully derived. Consideration of the third order 
of the interaction range is the key factor that allowed pre- 
diction of a novel type of soliton solution. It also helped 
us to find new acceptable ranges of physical values, where 
the existence of soliton is possible for the boson branch of 
the solution in the boson-fermion mixture, compared to the 
system comprised of bosons only. 
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